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ABSTRACT 

The transition of a standard thin disk to a radiatively inefficient accretion flow (RIAF) is expected to occur, 
when its accretion rate m is lower than the critical value merit ('w = M/MEAd)- The RIAF is very hot, and it radiates 
mostly in the hard X-ray band (> 100 keV). Assuming that the accretion disk in every bright active galactic 
nucleus (AGN) will finally undergo a RIAF phase while m < liicrh, we calculate the contribution of the RIAFs 
in AGNs to the cosmological X-ray background of 10-1000 keV. We find that the timescale f^''^^ of the RIAF 
accreting at ~ rhcnt should be shorter than ~ lQ~^t^ if mait = 0.01, where is the lifetime of bright AGNs, i.e., m 
declines from wicrit to a rate significantly lower than merit within t^^. The derived timescale t^^^^ is affected by 
the parameters adopted in the model calculations, which is also discussed in this Letter. 



Subject headings: galaxies: active — quasars: general- 
diffuse background 



-accretion, accretion disks — black hole physics; X-rays: 
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1. INTRODUCTION 

The AGN X-ray luminosity function (XLF) is directly linked 
to the accretion history of AGNs in the universe. Many works 
on XLFs were carried out in either the soft X-ray band (< 
3 keV) (e.g., Maccacaro et al. 1991; Boyle et al. 1993; Page et 
al. 1997; Miyaji et al. 2000) or the hard X-ray band (> 2 keV) 
(e.g., Boyle et al. 1998; Cowie et al. 2003; Ueda et al. 2003). 
The luminosity functions (LFs) derived from the surveys in the 
soft X-ray band (< 3 keV) may have missed many obscured 
(type II) AGNs, while the hard X-ray surveys (^ 2-10 keV) 
can trace the whole AGN population including obscured type II 
AGNs. The cosmological X-ray background (XRB) is mostly 
contributed by AGNs (Hasinger 1998; Schmidt et al. 1998). In 
the most popular synthesis models of the XRB based on the 
unification schemes for AGNs, the cosmological XRB from 
< 2 keV to more than several hundred ke V can be fairly well re- 
produced by using a template spectrum of AGNs consisting of 
a power-law X-ray spectrum with an exponential cutoff around 
several hundred keV (e.g.. Matt et al. 1994; Madau 1994; Co- 
mastri et al. 1995; Gifli et al. 1999; Ueda et al. 2003). Di 
Matteo & Fabian (1997) alternatively proposed that the very 
hard X-ray background (VHXRB, the term is used here for 
the hard X-ray background above 10 keV to distinguish from 
the conventionally mentioned hard X-ray background in 2-10 
keV) may be dominated by the thermal bremsstrahlung emis- 
sion from the advection dominated accretion flows (ADAFs) in 
low-luminosity AGNs. Further detailed ADAF spectral calcu- 
lations (Di Matteo et al. 1999) showed that many sources at 
redshiftz ^2-3 with ADAFs accreting at the rates close to the 
critical value are required to reproduce the observed VHXRB 
spectral shape. Recently, Ueda et al. (2003) derived a hard X- 
ray luminosity function (HXLF) from a highly complete AGN 
sample (2-10 keV), which includes both type I and type II 
AGNs (except Compton thick AGNs). Based on this HXLF, 
their synthesis models can explain most of the observed XRB 
from the soft X-ray band to the hard X-ray band around several 
hundred keV. Their calculations slightly («10-20%) underes- 
timate the relative shape of the XRB spectrum around its peak 
intensity (see Fig. 18 in their work). Such a discrepancy can be 



explained provided the same number of Compton-thick AGNs 
with logA^H = 24-25 as those with logA^H = 23-24 is included. 
This implies that the contribution to the VHXRB from the ra- 
diatively inefficient accretion flows (RIAFs) in AGNs may be 
important, but not dominant, even if the contribution from the 
Compton-thick AGNs is not considered. 

There are a variety of studies exploring the evolution of 
AGNs based on either optical quasar LFs or XLFs, or both 
(e.g., Haehnelt & Rees 1993; Haiman & Menou 2000; Kauff- 
mann & Haehnelt 2000; Yu & Tremaine 2002; Wyithe & Loeb 
2003; Marconi et al. 2004). A common conclusion is that the 
timescale of AGN activities is short compared with the Hubble 
timescale, though the quantitative results on the bright quasar 
lifetime vary from ^ 10^ to ~ 10^ years for different investiga- 
tions. The standard optically thick accretion disks are present 
in bright AGNs, provided the accretion rate is high. The AGN 
activity may be switched off while the gases near the black hole 
are exhausted (see Narayan 2002, for a recent review, and ref- 
erences therein). While the accretion rate m (m = M/MEdd) 
declines below a critical value merit, the standard disk transits 
to a RIAF (e.g., Narayan & Yi 1995). The RIAF is opticaUy 
thin, very hot, and its spectrum is peaked at around several hun- 
dred keV. There are numerous observational evidences indicat- 
ing that the RIAFs are indeed present in many low-luminosity 
AGNs and in our galactic center Sgr A* (e.g., Narayan & Yi 
1995; Lasota et al. 1996; Gammie et al. 1999; Yuan & Nai-ayan 
2004). 

In this Letter, we will explore how the inefficient accretion 
history of AGNs is constrained by the VHXRB. The cosmolog- 
ical parameters fijvi = 0.3, f^A = 0.7, and Hq = 70 km s"' Mpc~^ 
have been adopted in this work. 

2. MODEL 

The HXLF given by Ueda et al. (2003) is so far most suit- 
able for our present investigation, as it includes both type I and 
type II AGNs (except Compton thick AGNs). Assuming the 
accretion disk in every bright AGN to transit to a RIAF while 
its accretion rate is low, we derive the co-moving space number 
density of these faint AGNs containing RIAFs from the HXLF, 
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provided the accretion rate evolution is known. Based on theo- 
retical spectral calculations for the RIAFs, strict constraints on 
the accretion history of AGNs can be achieved from the com- 
parison with the cosmological VHXRB in 10- 1000 keV. 

2. 1 . Space density of faint AGNs 

We assume that all AGNs described by the HXLF have stan- 
dard accretion disks (we will come to justify this assumption in 
the Discussion section). Hereafter, we refer to these AGNs with 
standard disks as "bright AGNs", while those AGNs with RI- 
AFs as "faint AGNs". The X-ray luminosity in 2- 10 keV can 
be converted to the bolometric luminosity by using an empir- 
ical relation: Lboi = fcmLx, where /cor = 100 is adopted (Elvis 
et al. 2002; Menci et al. 2004). The black hole mass density 
for bright AGNs in the co-moving space at redshift z can be 
calculated by 
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where $(Lx,z) is the HXLF given by Ueda et al. (2003), wi"™"' 
is the average dimensionless accretion rate for bright AGNs de- 
scribed by this HXLF, and Z-Edd,© = 1.38 x lO^^ergs s"^ is the 
Eddington luminosity for a black hole with solar mass. 

In this Letter, we assume that the black hole mass does not 
change significantly after the accretion mode transition, which 
is satisfied only if we consider the time after the transition 
is shorter than the Salpeter timescale, because the accretion 
rates of these faint AGNs are very low (ni < 10"^). The total 
monochromatic X-ray luminosity of all faint AGNs in unit of 
co-moving volume is given by 
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Lx(E)dLx(E), 



(2) 



where n^(Mbh,z) is the black hole mass function for faint AGNs, 
L^'(Mbh,£') is the X-ray luminosity of the RIAF accreting at 
in = Went- For simplicity, we employ a conventionally adopted 
assumption of a fixed lifetime t^ and the same light curve for 
all faint AGNs. For the RIAFs in AGNs, their spectra LxiE) 
depend almost Unearly on the black hole mass Mbh, and Eq. (2) 
can be re-written as 

l}i'\E,z) = j, J ^^-«f(Mbh,z)dMbh J lxiE,t)dt 
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(3) 



where lx(E,t) is the faint AGN light curve for a lO^Mo black 
hole, and it can be calculated provided the time-dependent ac- 
cretion rate mit) is known. The total number density N^iz) of 
faint AGNs is: 7^(z) = N^iz^/t^, where is the Uftime of 
bright AGNs, and the bright AGN number density N^(z) can 
be calculated from the HXLF. The black hole mass density 
Pbh(z) = A^(z)M^™', where Mg™"" is the average black hole mass. 
The average black mass for faint AGNs should be larger than 
that for bright AGNs, as the black holes in bright AGNs are still 



growing through accretion (a rough estimate can give the aver- 
age hole hole mass in faint AGNs being about twice of that in 
bright AGNs, if m^""' = 1 and ~ 10^ years). This leads to 



f h 

Pbhiz) = -uPbb(Z)fbh, 



(4) 



where /bh > 1 is the ratio of average black hole masses of faint 
AGNs to bright AGNs. The black hole mass density for faint 
AGNs can be calculated from the HXLF by using Eqs. (I) and 
(4), so that Eq. (3) can be re-written as 
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Based on the RIAF models, the X-ray light curve can be calcu- 
lated, if we know how the accretion rate m evolves with time. 
Unfortunately, we are still ignorant of the detailed form of m(t). 
Here, we introduce a timescale t^^ to describe the basic fea- 
ture of the evolution of RIAFs in these faint AGNs, 
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This timescale f**^ describes how fast the accretion rate of a 
RIAF declines from merit to a rate significantly lower than rhcnt 
after the accretion mode transition. Now, Eq. (5) becomes 
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The contribution from the RIAFs in all faint AGNs to the 
cosmological XRB can be calculated by 



/- Lf [(l+z)£,z] dV^ 
J 4'Kdlil + z) dz 



1 



/bh rpl^(z)lf[(l+z)E]dV 



dz. (8) 



IO^Mq J 4ndlil+z) dz 



In this Letter, we conservatively adopt /bh = L 

2.2. RIAF spectra 

In order to calculate the contribution of RIAFs in faint AGNs 
to XRB, we need to have the X-ray spectrum lx\E) of a RIAF 
around a lO^^M^i black hole accreting at the critical rate merit. 
We employ the approach suggested by Manmoto (2000) to cal- 
culate the global structure of an accretion flow surrounding 
a Schwarzschild black hole (i.e., a = 0) in general relativistic 
frame. All the radiation processes are included in the calcula- 
tions of the global accretion flow structure (see Manmoto 2000, 
for the details and the references therein). Unlike Manmoto 
(2000)'s calculations, which are limited to the cases without 
winds, we also calculate the cases with winds. In the spec- 
tral calculations, the gravitational redshift effect is considered, 
while the relativistic optics near the black hole is neglected. 
This will not affect our final results on the XRB, as the faint 
AGNs should have randomly distributed orientations and the 
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stacked spectra would not be affected by the relativistic optics. 
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Fig. 1 . — The spectra of RIAFs accreting at the critical rate rhout = 0.01, for 
different wind parameters: pw = O(black), 0.2(green), 0.5(red), and 0.9(blue), 
respectively. The black hole mass Mth = IO^Mq, viscosity parameter a = 
0.2, the fraction of the magnetic pressure 1 — /3 = 0.5, and the outer radius 
Rout = 100/?schw are adopted in the calculations. The dotted lines represent the 
bremsstrahlung spectra of the RIAFs. 



Fig. 2. — The contribution to XRB from bright and faint AGNs for the 
RIAFs with different wind strengths or without winds. The dashed lines are 
the observed XRB. The lines with different colors represent different values 
of fl^'AF/^b^Q oo5(green), O.Ol(red), and 0.05(blue), respectively. The dot- 
dashed line represents the contribution by bright type I/II AGNs (Compton- 
thin), which is taken from Ueda et al. (2003). The solid lines represent the 
XRB contributed by bright type I/II (Compton-thin) AGNs and the RIAFs in 
the faint AGNs, while the dotted lines are for the contributions of RIAFs in the 
faint AGNs only. 



3. RESULTS 

The three-dimensional MHD simulations suggest that the 
viscosity parameter a in the accretion flows is ^0.1 (Ar- 
mitage 1998), or - 0.05-0.2 (Hawley & Balbus 2002). We 
assume a r-dependent accretion rate m = »iout('"/'"out)''* for the 
RIAFs with winds. The parameters adopted in the calcula- 
tions are: a = 0.2, ;7iout = 0.01, the fraction of magnetic pres- 
sure 1 -/? = 0.5 {f3 = Pgas/ptot), and the fraction of dissipated 
energy directly heating electrons is: S = 0.5. The outer ra- 
dius Tout = 100/?schw is adopted in all our calculations, where 
^Schw = 2GMbh/c^. We plot the X-ray spectra for RIAFs in Fig. 
1. 

There is no doubt that the contribution from a normal 
bright AGN is important in the VHXRB, as the BeppoSAX 
observations showed that the power-law X-ray spectra of 
bright AGNs extends to several hundred keV (e.g., Nicas- 
tro et al. 2000). Here we consider that the XRB con- 
sists not only of the emission from type I/II bright AGNs 
(Compton-thin) described by the HXLF, but also of the emis- 
sion from RIAFs in those faint AGNs which are not in- 
cluded in the HXLF. In Fig. 2, we plot the results for 
the RIAFs with different wind parameters p^ for differ- 
ent ratios t^^^^ /t"" (/«'"'" = 0.1 is adopted for bright AGNs). 



4. DISCUSSION 

We find that all spectra have an energy peak at around sev- 
eral hundred to 1000 keV (see Fig. 1). Compared with those 
obtained by Di Matteo et al. (1999), our spectra have higher 
energy peaks. The reason is that we adopt (5 = 0.5, larger than 
theirs, which leads to higher electron temperatures of the ac- 
cretion flows. From Fig. 2, we find that t^^^^ < 0.05?'' is 
required from the comparisons of our theoretical calculations 
with the observed XRB, for any RIAFs either without winds 
or with strong winds. For the RIAFs without winds (/?„ = 0) 
or weak winds (/?„ = 0.2), we find that the hard X-ray emission 
alone from the RIAFs in faint AGNs have already surpassed the 
observed XRB, provided f'^'^^ = 0.05f''. If the hard X-ray emis- 
sions from both the bright type I/II AGNs and RIAFs in faint 
AGNs are considered, we find more strict constraints on the 
RIAF timescale: t^^^^ < 0.0k'' required for most cases (it be- 
comes t^^'^^ < 0.005?'' for the RIAFs without winds). Here, we 
have neglected the contribution from the Compton-thick AGNs. 
If the contribution from a Compton-thick AGN is included, the 
present derived RIAF accretion timescales will become even 
lower. 

Our present calculations are based on the assumption that 
all AGNs described by the Ueda's HXLF are bright, i.e., the 
standard bright accretion disks are responsible for their energy 
sources. Our RIAF spectral calculations show that L^ '"'*^'^^ = 
1.48 X lO'^'ergs s"' for p„ = 0, and 5.40 x lO'^^ergs s"' for 
Pv, = 0.9, respectively, if Mbh = lO^^M© and nioui = 0.01. As 
the lower luminosity limit of the HXLF extends to lO'^'^ergs 
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s~\ it means that only the faint AGNs accreting at the criti- 
cal rate with black hole masses > 2 x lO^M© may have ap- 
peared in this HXLF. The standard disks accreting at m > merit 
around a black hole with > 2 x IO^Mq have L^^'^^"'^ > 10^^ ^ 
ergs s"'. This implies that some RIAF counterparts of the 
bright AGNs with L|"^°'^''^ > 10"*^^ ergs s~' may be included 
in this HXLF. We can roughly estimate that the number ratio of 
these RIAF counterparts to bright AGNs with L^-iokeV > ^q42.5 
ergs s"' is: < t^^^^ /t^. Simply integrating the HXLF, we find 
that the ratio of the sources with L^""""'^ = 10"*^'^^**ergs s~' to 
aU AGNs described by this HXLF is about 0.13, which im- 
plies that less than a fraction ~ 0.13?''^/?'' of all sources 
with L^"""'''^ = 10"*''^~"***ergs s"' may have RlAFs. Therefore, 
the contribution to the XRB from those RIAF sources with 
^2-iokeV y lo^i^ergs s"' may be over-calculated, but at a very 
low level of < 0.1 3?'^'^^/?'', which is negligible and will not 
affect our main conclusions. 

We have assumed a constant lifetime and (t^^) for all 
bright and faint AGNs, which may not be true, as suggested by 
Hopkins et al. (2005). However, our results on the XRB only 
depends on the ratio t^^^^ /t^, so our conclusions will not be al- 
tered, even if each individual sources have different lifetimes, 
provided they have a similar time-dependent form of m(t). The 
resulted XRB from Eq. (8) depends on the value of rh'^''^'^ as 
oc l/m'"^'=' (see Eq. 1). McLure & Dunlop (2004) estimated 
that the average accretion rate m^^^'^ varies from 0.1 at z ^ 0.2 
to 0.4 at z ~ 2 from a large sample of SDSS quasars. If we 
adopt a larger m^"^ = 0.4, the derived RIAF timescales 
will be four times of the present values. The derived timescales 
are proportional to the bolometric luminosity correction factor 
/cor, and the uncertainty of /cor should not be very large, which 



will not affect our main conclusions. The emissions in hard 
X-ray bands from RIAFs are dominated by the bremsstrahlung 
emission and Comptonization of the bremsstrahlung photons 
(see Fig. 1), which are almost independent of /3 (magnetic field 
strength) adopted in the calculations. 

Our results indicate that the accretion rate must drop to a 
very low rate from the critical accretion rate in a short timescale 
compared with the bright AGN Ufetime. Recent hydrodynam- 
ical simulations on the evolution for the AGN formed after a 
merger indeed show a rapid decline of accretion rate from near 
Eddington rate to ~ 10"^ within a short timescale compared 
with its bright phase (see Fig. 2 in Di Matteo et al. 2005), 
which is qualitatively consistent with our results. The black 
hole growth during the timescale t^^'^^ can be neglected com- 
pared with that during the bright AGN phase, as Jriaf <^ t^- 
However, our present calculations have not considered how the 
accretion rate evolves with time in detail. The advection be- 
comes important while the RIAFs are accreting at rates <C Wcrit, 
the black holes may swallow gases without radiating much X- 
rays. So, it cannot be ruled out the possibility that the faint 
AGNs stay at accretion rates far lower than the critical value for 
a very long time, say, comparable with the Hubble timescale. 
In principle, this can also be constrained by the cosmological 
XRB, which is beyond the scope of this Letter. 
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